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SUMMARY 


Flight tests with the Transport Systems Research Vehicle (TSRV) Boeing 737 air- 
plane have shown that advanced forms of guidance presented on electronic flight dis- 
plays may result in improved lateral and vertical tracking during the approach to 
landing. A critical element in producing the advanced guidance is the accuracy of 
the position estimate computed from inputs from radio navigation signals in the TSRV 
area navigation system. A comparison was made of the accuracy of two position esti- 
mation algorithms, each using a combination of instrument landing system (ILS) local- 
izer and distance measuring equipment (DME) radio signals from an arbitrarily located 
DME facility to compute new estimates. One algorithm computes a position estimate 
with two independent error components, one from DME and the other from ILS localizer 
signals. The other algorithm computes a position estimate with a single error com- 
ponent, which is found by geometrically combining the output of both DME and ILS 
localizer signals. 

Both simulation tests and flight tests using a radar tracker to determine flight 
path accuracy showed that an artificially induced position estimate error was reduced 
faster with the single-component solution than with the dual-component solution. 
Faster time constant gains were used in the position estimate filter for both algo- 
rithms to reduce position errors faster. A NASA test pilot reported that flight 
guidance based on the position estimate was smooth for all algorithm configurations 
and filter time constants tested. 

A bias error was noted in the longitudinal axis (relative to the runway). this 
error was caused by a bias error within the DME indication. However, the bias error 
was well within the tolerances specified for the en route DME systems used during 
these tests. Hie resulting error may require the use of precision DME for guidance 
during the approach to landing segments of flight. 


INTRODUCTION 

Flight tests with the Transport Systems Research Vehicle (TSRV) Boeing 737 air- 
plane flown in the microwave landing system (MLS) environment have shown that the 
TSRV area navigation system can provide situation and guidance information, in the 
form of a cathode ray tube (CRT) map display and a computer-drawn perspective runway 
symbol on the primary flight display (fig. 1), which results in improved pilot lat- 
eral and vertical tracking during the landing approach (refs. 1 and 2). A critical 
element in producing this advanced guidance was the accurate estimate of airplane 
position computed from the MLS signals by the navigation system. However, MLS is not 
yet commonly available at most airports; thus other forms of radio navigation aids 
will be required to determine airplane position if an earlier implementation of these 
advanced display concepts is desired. 

The instrument landing system (ILS) is a precision radio navigation aid avail- 
able at many airports. However, the ILS provides lateral information in the hori- 
zontal plane in the form of angular deviation relative to the runway centerline. 

This information must be combined with range information to completely define air- 
plane horizontal position relative to the runway. An initial attempt was made to 
compute the position estimate with an algorithm that determined two independent 



position difference components, one from ILS localizer signals and the other from DME 
signals. The ILS-based component was perpendicular to the runway centerline and the 
DME-based component was on a line between the airplane and the DME ground station. 
These components were used as inputs to the position estimate filter in the area 
navigation system. Although the design used in this first algorithm resulted in a 
more accurate position estimate than was normally obtained with the dual-DME position 
estimate mode (ref. 3), it did not provide sufficient accuracy for the pilots to use 
the advanced display symbology for primary guidance during an approach to landing. 

To improve the accuracy of the position estimate, a second algorithm (originally 
documented in ref. 4) was designed for the TSRV 737 airplane navigation system. This 
algorithm used both ILS and DME information to compute a single position difference 
component for the position estimate filter. The locations of the ILS localizer and 
DME antennas were stored in the navigation data base and azimuthal information from 
the ILS localizer was geometrically combined with distance information from an arbi- 
trarily located DME station to calculate position. The output of this geometric 
calculation was used as an input to the position estimate filter to produce a 
smoothed estimate on which display guidance could be based. 

Simulation tests were conducted to adjust the filter time constants so that 
large position estimate errors occurring during en route flight would be nulled 
quickly while maintaining a smooth display of guidance to the pilot during transition 
into the ILS environment. Flight tests were conducted with the TSRV 737 airplane in 
a ground-based radar tracking environment to document the difference between the 
actual position of the airplane (radar defined) and the area navigation system posi- 
tion estimate. The purposes of this report are to document the equations and logic 
used to generate the position estimate with the ILS and DME information and to dis- 
cuss the results of these simulation and flight tests. 


SYMBOLS 

A distance between the DME and ILS localizer antenna, n.mi. 

a east component of A, n.mi. 

b north component of A, n.mi. 

D t distance from DME to airplane corrected for slant range, n.mi. 

DJ_ slant range distance from DME to airplane, n.mi. 

DP position difference vector 

DP component of position difference vector DP perpendicular to the runway 

" centerline 

DPp magnitude of vector DP^, n.mi. 

F ellipticity constant, 0.003367 

h altitude of airplane above mean sea level, ft 

cl 

h nM „ altitude of DME antenna above mean sea level, ft 

UMd 
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X' , 
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altitude of ILS localizer antenna above mean sea level, ft 
unit coordinate vectors 
position estimate filter gain 

_ <1 

velocity estimate filter gain, hr 
velocity estimate filter gain 

outer radial limit of ILS localizer geometric check, n.mi. 

inner radial limit of ILS localizer geometric check, n.mi. 

axes of orthogonal coordinate system oriented toward true north 

angle formed by the vector Z and a line between the airplane and DME 
antenna, deg 

radius of Earth, ft 

meridional radius of curvature, n.mi. 

normal radius of curvature, n.mi. 

unit vector perpendicular to runway centerline 

north, east components of navigation system velocity estimate, 
knots 

north, east components of inertial ground speed, knots 

axes of orthogonal coordinate system oriented along the runway centerline 
coordinates of vector Z^ transformed into the X’,Y' coordinate system 
vector of estimated airplane position from ILS localizer antenna 


e,N' 


m 


m 


*m,N' 


magnitude of Z , n.mi. 

e 

Z _ north, east components of Z , n.mi. 

e / e 

vector of measured position of airplane from the ILS localizer antenna 

->■ 

magnitude of Z , n.mi. 

m 

Z „ north, east conponents of position vector Z , n.mi. 

in, r. m 


Z vector of estimated airplane position measured radially along ILS deviation 

angle 

Z_ magnitude of position vector Z , n.mi. 

r r 

Z r,N' Z r,E north, east components of position vector Z , n.mi. 
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a angle formed at ILS localizer antenna by DME antenna and measured airplane 

position, deg 

P angle P limit at which DME inputs will be inhibited, deg 

AP ,AP north, east components of position estimate difference, n.mi. 

N E 


At 


computer iteration cycle, sec 


AV n ,AV e north, east components of system velocity update, knots 
Acp , AX latitude, longitude update terms, deg 

r) localizer deviation angle relative to the runway centerline, deg 

p relative angle between the DME and the ILS localizer antenna, deg 

latitude, longitude of DME antenna location, deg 
latitude, longitude of estimated airplane position, deg 
latitude, longitude of ILS localizer antenna, deg 


^DME' \>ME 

<t> » ^ 

e e 


t> 

o o 




Q 


angle of runway centerline (opposite to direction of landing) relative to 
true north, deg 


vertical angular limit of ILS localizer geometric check, deg 
a) lateral angular limit of ILS localizer geometric check, deg 

Subscripts: 
t at time t 

t-1 at iteration time previous to time t 

Abbreviations: 

ACY Atlantic City DME navigation radio facility 

cathode ray tube 
distance measuring equipment 

navigation position estimate mode: inertial velocity and dual DME 


CRT 

DME 

IDD 

ILD 


navigation position estimate mode: inertial velocity, ILS localizer, 

and DME 

ILS instrument landing system 

I LX navigation position estimate mode: inertial velocity and ILS localizer 

INS inertial navigation system 
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MSL microwave landing system 

NCU navigation computer unit 

rms root mean square 

TSRV Transport Systems Research Vehicle 

VCN Cedar Lake DME navigation facility 

VOR very-high-frequency omnidirectional-range navigation radio 


DESCRIPTION OF NAVIGATION POSITION ESTIMATE 
General Description 

The TSRV navigation computer was programmed to compute an estimate of airplane 
position and velocity based on combining data from an inertial navigation system 
(INS) and various radio navigation systems. Information from the airborne central 
air data computer and the magnetic direction system is also available in the event 
that INS information is unavailable. 

Signals from the radio navigation systems are combined to determine airplane 
position. However, because of signal noise and various errors inherent within the 
radio navigation systems, this definition of airplane position must be combined with 
inertial velocity information from the INS in a position estimate filter to produce a 
smoothed estimation of airplane position. This estimate of position is computed at 
20 Hz and is used as a basis for navigation, guidance, and display computations. 

Various position estimate modes may be used, depending upon availability of the 
radio navigation signals. Position estimates normally are calculated with inputs 
from two different DME facilities for flights away from the terminal area. This 
position estimate mode is identified by the letters IDD (inertial velocity, DME, DME) 
shown on the electronic map display in front of each pilot. While the aircraft is 
operating in the terminal area and within signal coverage of an instrument landing 
system (ILS), the position estimate will be computed either with both ILS localizer 
and DME radio signals (ILD) or, if DME is not available, with localizer signals only 
( ILX) . 


In the TSRV navigation system, the radio signals are used to conpute a position 
difference term to be used in the position estimate process. The position difference 
term is defined as the vector between the last previously estimated position and the 
position defined by the radio navigation information. The position difference vector 
is divided into north and east components, which are used as inputs in the position 
estimate filter. The equations used in the position estimate filter are given in 
appendix A. 


ILS-DME (ILD) Position Estimate: Dual Component 

The design of the first ILS position estimate algorithm involved simply substi- 
tuting an ILS localizer input for one of the DME inputs in the conventional dual-DME 
position estimate algorithm. This resulted in two independent position difference 
components. The equations for the difference component computed with localizer 
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information are also used in the ILX mode and are described in appendix B. Hie equa- 
tions for the difference component computed with the DME signal are described in 
reference 3. Each of these difference terms is resolved into north and east compo- 
nents, summed, and then used as input in the position estimate process. 


ILS-DME (ILD) Position Estimate: Single Component 

General solution. - Figure 2 shows the geometry of the ILS localizer antenna, the 
DME antenna, the airplane estimated position, the airplane position measured with 
localizer and DME information, and the runway with an extended centerline. An 
orthogonal coordinate system, with its origin placed on the localizer antenna, is 
oriented with respect to true north. Since the DME could be selected either manually 
or automatically, the position difference equations were derived so that the DME was 
not required to be collocated with the ILS localizer antenna. 

By determining the relative geometry between the ILS localizer antenna, the DME 
antenna, and the position of the airplane measured with both localizer and DME infor- 
mation, the vector Z between the origin and the measured position of the 

m, t 

-> 

airplane was found. Hie vector Z „ was found directly from the known latitudes 

e , t-1 

and longitudes of the origin and the previously estimated position of the airplane. 

-> ->• 

The positron difference was then found by subtracting Z from Z and 

m,t e,t-1 

dividing into the north and east components AP n and AP^ 

-y -y 

Calculation of airplane position vector Z .- Figure 2 shows that Z is 

m, t m, t 

one side of the triangle formed by the origin (ILS localizer antenna), the measured 
position of the airplane, and the DME antenna. Known quantities used to determine 

Z in this triangle include the angle of the runway centerline (opposite to 

m, t r 

I 

the direction of landing), the slant range D t measured from the DME signals, the 

localizer deviation angle T) t relative to the runway centerline and measured from 

the ILS localizer signals, the location (latitude and longitude) ^dme'^DME of the 

DME antenna, and the location <}> ,X of the ILS localizer antenna. 

o o 

Figure 3 shows the angular geometry and distance A between the origin and the 
DME antenna. Hie length A and the angle |i between the DME antenna and the north 
axis remained constant for a particular DME and were calculated only once. If a 
different DME station was tuned, the A and \i were recalculated. Hie length of 
A was determined by vectorially summing its components a and b: 


a = 60(\ 


DME 


b = 60(6 


DME 


X ) cos 
o 


DME 


n.mi . 


6 ) n.mi. 
o 


6 



n.mi 


A 



+ b 


2 


The angle p was found by 


V- = 



(0 < p < 2n) 


The angle a formed by side A and vector may vary continuously and 

was calculated 20 r times per second. It ranged betweefi 0 and n, inclusive, and was 
found by 


\ = 1 2 it - p + (4> r - T) t )| 


If this calculated value of 


was greater than 


n, then 


= 1 2n - [\2% - p + <4> r - T l t )|]| 

I 

The value of the DME reading D fc measured in the airplane was the slant range 
distance between the ground-based DME antenna and the airplane. This distance was 
corrected to determine the ground distance D fc shown in figure 3 between the air- 
plane and DME: 


D 


t 



sin 


- 

/ h - h \ 


„ -1/ 

a, t DME' 


cos 

V D t ) 


_ 


- 


n.mi. 


Angle formed by the vector 

was calculated 20 times per second. 


Z and side D fc 

Frofti the relation 


could vary continuously and 


A 

sin P 


D 

t 


sin a 


n.mi. 


angle P fc was determined by 


P 


t 


. -1 
= sm 



A comparison of the square of side A with^the sum of the square of side D fc and 
the square of the magnitude of the vector Z was used to determine whether 
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angle P t was obtuse or acute. (The airplane position vector estimated on the last 

iteration (Z ) was used as an approximation for Z . .) Hence, 

e,t-1 “i 1 


r 




sin "Yf t sin a t) 


. -1/A 

n - sin f— — sin a 
t 


L 2 <o 2 + z t 2 ) 

\ t e,t-1 f 


) < 


2 2 
A > D + Z 
t e 


2N > 

,t-1 ) 


~r 

The magnitude of the estimated airplane position vector Z was 

0 | L“ I 


>,t-1 ^ ^ Z e,N,t-1 ^ + (Z e,E,t-1 ) ^ 


2 1 1 /2 


n.mi. 


in which 


Z = 60(4) - 4> ) n.mi. 

e,N,t-1 e,t-1 o 


Z .-I . . — 60( X ^ . 
e,E,t-1 e,t-1 o 


<b + 4> 

, , o e,t-1 

- X ) cos — 2 ~~ — n.mi. 


The magnitude of Z was then found as follows: 

m, t 


Z = A cos a + D cos P. n.mi. 
m,t t t t 


The north and east component^ of Z^ ^ were calculated from the angle <J> - 

between the north axis and Z : ' 

m, t 


Z = Z cos((l) - Ti ) n.mi. 

m,N,t m, t r t 


Z = Z sin((|) - H ) n.mi. 

m,E,t m,t r t 


Calculation of position difference in north and east components AP ,AP 

N j t E, t 

The position difference in north and east components was then found by subtracting 
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and subtracting 


-V 

the north component of z 

e,t-1 

-y 

the east component of Z 

6 f i 


from the north component of Z 

m, t 

from the east component of Z : 

m, t 


AP 


N, t 


= Z 


m, N, t 


- Z 


e,N, t-1 


n.mi. 


AP 


E, t 


= Z 


m, E, t 


- Z 


e,E,t-1 


n.mi 


These position difference components were then used directly in the navigation posi- 
tion estimate algorithms previously described. 


ILS Localizer-Airplane Position Validity Check 

Since random DME selection, automatic frequency tuning, and other means of auto- 
matic software control were utilized in the navigation computer, a check had to be 
made to ensure that the navigation computer was using the appropriate navigation data 
for the area in which the airplane was flying. Obviously, if improper navigation 
data were being utilized, position estimates based on those data would be inaccurate 
and should not be used. 

To preclude the use of the wrong ILS localizer navigation data, a geometric 
check was made to determine whether the estimated position of the airplane was within 
the localizer boundaries. If the airplane was not within these boundaries, then ILS 
updating was inhibited. These boundaries included lateral and radial limits of 
coverage, as shown in figure 4. The lateral limit of the localizer angle of coverage 
is +oj from the runway centerline. Radial limits required that the airplane be 
within distance L of the localizer antenna, but not closer than distance M. A 
vertical angle of coverage limit Q was measured from the localizer antenna. 

A new orthogonal coordinate system, with its origin located at the ILS localizer 
antenna and its X'-axis on the runway centerline, was used to make the geometric 
check. The north and east components of the estimated position vector Z were 

transformed into the new X' ,Y' coordinate system as follows: e ' 


X' 

= Z 

COS 

<i> 

+ Z 

e. 

t e,N, t-1 


r 


Y> 

= z 

sin 

<\> 

- z 

e. 

t e,N, t-1 


v r 



The lateral azimuth check was then 

-X' ^ tan 0 ) < Y' < X* 
e,t e,t e,t 


e, E, t-1 


e, E, t-1 


tan oj 


sin (); 

r 

cos cl 

r 

n.mi . 


n.mi. 


n.mi . 
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The radial check was 


X' 

0 t 

M < L < L n. mi. 

cos n t 


and the vertical check was 


h < h + X' tan Q 
a,t o e, t 

For these tests, the ILS localizer azimuth coverage limits were 


Lateral, to, deg 20 

Vertical, Q, deg 10 

Inner radial, M, n.mi 0.165 

Outer radial, L, n.mi 10 


A geometric check of the location of the DME ground station relative to the 
airplane and the runway centerline is performed continuously to prevent the possibil- 
ity of large position estimate errors caused by the DME inputs (fig. 5). As the 
angle P fc formed by the localizer antenna, airplane, and DME ground station 
approaches 90°, small errors from the DME range information can result in signifi- 
cantly larger radial position errors along the localizer deviation angle on which the 
airplane is located. Hence, if angle falls within a prescribed angular boundary 

(3 about P t = 90°, the DME signals will not be used and a new DME will be tuned or 
I LX updating will occur. This geometric check is defined as follows. If 


(f + p) > h > (I - <>) 

then do not use DME information. For the simulation and flight tests documented in 
this report, f3 = it/4 or 45°. 


DESCRIPTION OF SIMULATION, AIRBORNE FLIGHT SYSTEMS, 

AND GROUND-BASED RADAR 

TSRV 737 Simulation 

Hie TSRV 737 simulation consists of a fixed-base cockpit identical to the 
research flight deck in the TSRV airplane and a sophisticated, six-degree-of-f reedom 
digital representation of the airplane and its experimental systems. Nonlinear 
effects such as engine response, varying stability functions, and control servo 
models, as well as modeling of the navigation, guidance, and display subsystems, are 
used to represent the research airplane with a high degree of realism. 
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Airborne Flight Systems 


Airplane. - The NASA test airplane is the Boeing 737-100 twin-jet transport air- 
plane shewn in figure 6. Hie TSRV airplane is used as a research vehicle and has a 
complete set of experimental navigation, guidance, flight control, and display sys- 
tems in a separate research cockpit located in the cabin of the airplane. All the 
normal flight systems (e.g. , flight control, navigation, pressurization) have been 
retained in the conventional cockpit in a normal, functional state. This arrangement 
allows changes to be made to any of the experimental systems while retaining the 
standard operational features of the airplane. 

Experimental systems. - The ejqperimental research flight systems consist of a 
triplex digital flight control computer system, an electronic CRT display system, and 
a digital navigation and guidance system integrated into a separate, two-man-crew, 
research cockpit {ref. 3). The navigation and guidance calculations are performed in 
a single digital navigation computer. Various navigation sensor signals (including 
inertial navigation system velocities and accelerations, true airspeed, magnetic 
heading, and VOR, ILS, and DME radio signals) are used in the navigation computer. 
This computer estimates the airplane position based on combinations of the sensor 
inputs. Horizontal path, vertical path, and thrust guidance and commands based on 
the estimated position, estimated velocity, and path tracking errors of the airplane 
are displayed to the pilot on the CRT display system. 

Airbor ne data acquisition system .- A wideband magnetic-tape recorder onboard the 
airplane recorded data at a rate of 40 samples per second. These data included 
93 parameters describing the airplane configuration, attitude, and control surface 
activity. Thirty-two additional channels were used for recording navigation computer 
parameters, including the latitude and longitude of the navigation position estimate. 
Video recordings of the electronic attitude director and electronic horizontal- 
situation displays were also made throughout the flights. 


Radar Tracking Facilities 

Radar tracking of the airplane was provided by the Federal Aviation Administra- 
tion Extended Area Instrumentation Radar (EAIR) Facility at the Federal Aviation 
Administration Technical Center (FAATC) in Atlantic City, New Jersey. The tracking 
radar is a precision C-band instrumentation radar system that was operated in a 
beacon-tracking mode during these flight tests. Slant range, azimuth angle, and 
elevation angle data were recorded on magnetic tape at a sample rate of 10 Hz. All 
airborne and ground-based recorded data were time correlated for postflight process- 
ing and analysis. Postflight processing of the range, azimuth, and elevation data 
consisted of conversion to latitude, longitude, and altitude. 

Radar accuracy is stated (ref. 5) as 0.2 mrad, (0.011°) in azimuth and elevation 
angle and less than 20 yards rms error in range. 


RESULTS AND DISCUSSIONS 
Simulation Tests 

Test procedure. - The TSRV 737 airplane simulation was used both in a fast-time 
mode without the research cockpit and in a real-time mode with a pilot operating the 
research cockpit. The fast-time mode was used to investigate the change in the 
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position estimate convergence characteristics for the single- and dual-component 
position difference solutions and to select gains for the position estimate filter. 
The real-time piloted mode was used to determine the effects that airplane maneuvers 
might have on the guidance based on the navigation position estimate. 

Figure 7 is a plan view of the reference coordinate system, the runway, the DME 
antenna, and the localizer antenna modeled for the simulation. Hie origin of the 
reference coordinate system was located at the approach end of the runway with the 
x-axis oriented along the centerline of a north-south runway. The ILS localizer 
antenna and the DME antenna were collocated on the x-axis 8500 ft past the runway 
threshold (reference axis origin). This geometrical arrangement of the reference 
axes, runway alignment, and antenna locations was selected for these tests to facili- 
tate data analysis. With this geometrical arrangement, the lateral component (east- 
west) of the position estimate was derived from ILS localizer signals and the lon- 
gitudinal component (north-south) was derived from DME signals. 

Fast-time simul ation.- During the fast-time simulation tests, the airplane was 
constrained to fly along the ILS localizer toward the runway regardless of the navi- 
gation position estimate. At the start of each test, the navigation position esti- 
mate was offset to create a position estimate error of 1000 ft to the right and 
1000 ft to the rear of the actual position of the airplane. The subsequent reduction 
of these errors was used to document the convergence characteristics for each test 
run. Hie position estimate error was broken into lateral and longitudinal components 
relative to the localizer centerline so that convergence due to ILS localizer infor- 
mation (lateral) and DME information (longitudinal) could be analyzed separately. 

The convergence characteristics for each component were quantified by the initial 
convergence rate, the time required to reduce the initial displacement error by 
63 percent, and the resulting amount of overshoot. 

The first four fast-time simulation runs were conducted to compare the conver- 
gence characteristics of the single-component solution with those of the dual- 
component solution. Each of these solutions was run with position estimate filter 
gains based on the 50 sec time constant normally used with dual-DME estimates (test 
condition 1) and with a faster, 30-sec time constant (test condition 2). The result- 
ing lateral and longitudinal convergence characteristics for these runs and test con- 
ditions are tabulated in table I. 


TABLE I.- LATERAL AND LONGITUDINAL CONVERGENCE CHARACTERISTICS OF ILS-DME POSITION ESTIMATE 

[TSRV simulation; = 1 ] 





Initial recovery 

63-percent recovery 

Offset at 


Test 

Component 

time 

rate, 

fps 

time 

sec 

threshold,* ft 

Radio 

relative to 







signal 

condition 

constant. 


Single 

EUal 

Single 

Dual 

Single 

localizer 

sec 

Dual 

noise 



component 

component 

component 

component 

component 

component 


1 

Lateral 

50 

20.5 

38.0 

45.0 

23.7 

-61 .7 

-46.9 

No 

2 

Lateral 

30 

33.9 

64.3 

27.0 

14.2 

-80.8 

-35.7 

No 

1 

Longitudinal 

50 

40.1 

39.1 

22.8 

23.6 

-49.3 

-52.0 

No 

2 

Longitudinal 

30 

65.9 

64.4 

13.7 

14.1 

-45.6 

-42.9 

No 


♦Negative offset indicates an overshoot. 



Table I shows that little change occurred in the longitudinal convergence char- 
acteristics between the single- and dual -component position difference solutions. 

This was expected because of the collocation of the DME and localizer antennas. In 
this geometric configuration, the magnitude of the position difference term was the 
same for each solution. Hence, the convergence characteristics were the same. Use 
of the 30-sec gains resulted in approximately 66 percent faster convergence than was 
obtained with the 50-sec gains. 

The lateral convergence characteristics shown in table I indicate that use of 
the single-component position difference solution resulted in a faster reduction of 
error than was obtained with the dual-component solution. These results are illus- 
trated in figure 8, which shows a time history trace of the initial 1000-ft lateral 
error and the subsequent recovery for each of the four runs. Use of the single- 
component solution resulted in an 85-percent increase in the initial lateral recovery 
rate (90 percent with the 30-sec time constant) and a 48-percent decrease in time to 
reduce the lateral error by 63 percent of the initial offset (47 percent with the 
30-sec time constant) . Overshoot occurred on all four runs. The overshoot was 
reduced from 80.8 ft with the dual -component solution to 35.7 ft with the single- 
component solution (30-sec time constant). However, it was felt that the overshoot 
occurring on all runs could result in flight crew confusion due to disagreement 
between guidance based on the navigation position estimate and guidance based on 
localizer deviation information. 

The computations used to calculate the navigation system velocity estimate were 
examined to determine their effects on the position estimate convergence. The system 
velocity estimate is calculated by summing ground speed from the inertial navigation 
system with a comparatively small velocity correction term (AV N and Av e fc ) com- 
puted from the position difference defined by navigation radio information . ' This 
correction term is computed by summing a portion of the current position difference 
to the value that was computed on the last computation cycle (K 3 = 1.0) as follows: 


AV 

= K 

AV 

+ K 

ap 

knots 

N, t 

3 

N, t-1 

2 

N, t 


AV 

= 

AV 

+ K 

AP ^ 

knots 

E, t 

3 

E, t-1 

2 

E, t 



Bias errors from the INS velocity inputs may be nulled over a period of time by 
accumulating a portion of the position error caused by the velocity error. However, 
when a large position step error is suddenly induced, as might be obtained when 
changing from en route navigation radio facilities to approach facilities, the veloc- 
ity correction term will become great and could cause the navigation position esti- 
mate to momentarily undershoot or overshoot during the subsequent position error 
nulling process. The large position step error would eventually be nulled, but this 
would require more time than necessary to complete the approach to landing. In an 
effort to reduce the transient effects of large step errors the position estimate, 
the accumulation feature of the velocity estimate equation was removed by setting 

K 3 = °* 

The four fast-time simulation runs (test conditions 1 and 2) were repeated using 
the system velocity estimate equations with set equal to zero. The convergence 

characteristics resulting from these runs are quantified in table II. Figure ‘9 shows 



TABLE II.- LATERAL AND LONGITUDINAL CONVERGENCE CHARACTERISTICS OP ILS-DME POSITION 
ESTIMATE USING MODIFIED SYSTEM VELOCITY ESTIMATE EQUATIONS 

[TSRV simulation; K 3 = 0] 





Initial recovery 

63-percent recovery 

Offset at 


Test 

Component 

time 

rate 

fps 

time 

sec 

threshold,* ft 

Radio 

condition 

relative to 
localizer 

constant. 

Dual 

Single 

Dual 

Single 

Dual 

Single 

rid 1 

noise 




component 

component 

component 

component 

component 

component 


1 

Lateral 

30 

33.8 

64.1 

29.2 

14. 8 

3.2 

0.01 

No 

2 

Lateral 

50 

20.4 

39.0 

48.7 

24.7 

33.5 

1 .0 

No 

3 

Lateral 

30 

32.9 

62. 3 

31.1 

16.2 

13.7 

8.5 

Yes 

4 

Lateral 

50 

19.9 

37.9 

51 .4 

26.5 

44.5 

1 0.7 

Yes 

1 

Longitudinal 

30 

65.8 

64. 3 

14.3 

14.7 

8.7 

7.0 

No 

2 

Longitudinal 

50 

40.0 

39.1 

23.8 

24.6 

8.5 

6.9 

No 

3 

Longitudinal 

30 

65.9 

64.5 

14.2 

14.6 

4.4 

2.7 

Yes 

4 

Longitudinal 

50 

40.2 

39.2 

23.6 

24.4 

2.5 

0.3 

Yes 


♦Positive offset indicates an undershoot. 


the time history trace of the lateral component for each of these runs. These data 
indicate that the overshoot was substantially reduced. However, undershoot resulted 
in the run using the dual-component solution and the 50-sec filter gain. More time 
was required for nulling the lateral error with this solution and gain. The initial 
recovery rates were about the same as those obtained with the unmodified velocity 
estimate equations. The times required to reduce the initial displacement error by 
63 percent were increased slightly due to the better estimate of system velocity 
computed with the modified equations. 

Randomly generated noise was added to the simulated ILS localizer and DME radio 
signals to note any adverse effects on the convergence characteristics, particularly 
those obtained with the faster time constants. The four fast-time simulation runs 
were then repeated using the modified system velocity estimate equations (test con- 
ditions 3 and 4) . The convergence characteristics resulting from these runs did not 
differ significantly from those without noise. The results are quantified in 
table II. 

Real-time simulation.- A NASA test pilot was used to evaluate changes to the 
approach guidance resulting from the various combinations of position difference 
solutions and time constants. The pilot's task was to fly an ILS approach to the 
runway using the computer-drawn runway on the primary flight display (fig. 1 ) and 
associated symbology (described in ref. 2) for guidance. The pilot was to subjec- 
tively evaluate changes to the guidance resulting from modifications to the position 
difference solutions and the position estimate time constants. To aid the pilot with 
this evaluation, a second computer-drawn runway was drawn on the primary flight dis- 
play to represent the real-world runway as it would appear through a superimposed 
display of a forward-looking television image. The real-world runway had no simu- 
lated noise or position errors and was to be used for comparison with the computer- 
drawn runway based on the navigation position estimate. 

The evaluation runs were made with both single- and dual-component difference 
solutions, both 30- and 50-sec position estimate time constants, and both modified 
(K^ = 0) and unmodified (K^ = 1) velocity estimate equations. Simulated radio noise 
was also included during these runs. 




One set of runs was made with an initial lateral position estimate error of 
1000 ft. This position offset was readily seen in the difference between the real- 
world runway and the computer-drawn runway. The pilot did not report any difficul- 
ties or confusion in using the approach guidance while the position estimate 
algorithm was nulling the position error. However, both the overshoot, which 
resulted from the use of the unmodified velocity estimate equations, and the under- 
shoot, which was obtained with the dual-component solution, modified velocity equa- 
tions, and 50-sec gains, were very apparent, as could be seen by the differences in 
alignment between the simulated real-world runway and the computer-drawn runway. The 
pilot reported that the misalignment of the computer -drawn runway would degrade its 
usefulness when flying closer than approximately 1 mile to the runway. The use of 
the single -component difference solution, the modified velocity equations, and the 
30-sec time constants resulted in good runway alignment and provided adequate guid- 
ance information for flight to the runway threshold. 

An additional set of runs were made with no initial position estimate error. 

The pilot was instructed to turn the airplane back and forth approximately 20° to 
either side of the runway centerline to judge the dynamics of the display symbology 
caused by track angle and lateral displacements of the airplane. The pilot judged 
the motion of the symbology during these maneuvers to be smooth and acceptable for 
approach guidance for all runs. 


Flight Tests 

Te st procedure .- Flight tests were conducted to validate the results obtained 
from the simulation tests. The TSRV 737 airplane was flown on a series of ILS 
approaches to runway 1 3 at the Federal Aviation Administration Technical Center 
(FAATC ) airport located near Atlantic City, N.J. Ground-based precision radar track- 
ing data were recorded to define the actual position of the airplane during each test 
run. During postflight data analysis, the vectorial difference between the radar- 
defined position and the onboard estimated position was computed and resolved into 
longitudinal (parallel to the runway centerline) and lateral (perpendicular to the 
runway centerline) components. The locations of the DME's used during these tests 
(fig. 10) resulted in the lateral difference component being influenced primarily by 
the ILS localizer signals and the longitudinal difference component being influenced 
by the DME signals. 

Figure 10 shows the flight path used during these tests. The path was pro- 
grammed into the area navigation system to provide tracking guidance for the pilot to 
intercept the ILS localizer at a point 10 to 15 n.mi. from the localizer antenna. 

Position estimates computed during the flight to the final approach path used iner- 

tial velocity and either dual- or single-DME inputs, depending on the signal and 
geometrical validity requirements programmed into the navigation system. The algo- 
rithms used during this process are described in reference 3. 

Once the ILS localizer was intercepted, the pilot would slow the airplane to 
approximately 140 knots and maintain level flight at an altitude of 3500 ft. Radar 

and airborne data recording were begun prior to passing waypoint ILS10, which was 

located 10 n.mi. from the localizer antenna. After passing waypoint ILS1 0, the navi- 
gation position estimate was computed using the ILS and DME information. Once the 
glide slope was intercepted, the pilot used ILS localizer and glide slope deviation 
information for guidance since these were unaffected by the position estimate errors 
artificially induced for these flight tests. Each test run was complete once the 
airplane crossed the approach threshold of the runway. 
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During the test run, when the airplane was approximately 8 n.mi. from the local- 
izer antenna, the flight test engineer inserted an incremental change in the position 
estimate which resulted in a 1000-ft right lateral shift (relative to the ILS local- 
izer) and a 1000-ft rearward longitudinal shift. Hie convergence characteristics 
resulting from the correction of these induced errors were quantified with the same 
parameters used in the simulation tests. 

A subjective evaluation was made by a NASA test pilot to determine if the guid- 
ance based on the navigation position estimate was adversely affected by the modified 
estimate algorithms and the faster filter time constants. The evaluation was based 
primarily on a comparison of the computer -drawn runway with the real-world runway 
television image superimposed on the primary flight display. 

Functional check. - An initial test run was flown to conduct a functional check 
of the airborne software and to determine that adequate DME signal coverage was 
available for these flight tests. During this initial test run, ILS localizer and 
DME signals were used 93 percent of the time to compute the position estimate. 

During the remainder of the test run, inputs from the DME receiver were not valid and 
the position estimates were computed with localizer inputs only (ILX mode). 

The automatic tuning feature of the navigation system has been programmed to 
select DME radio stations that meet both geometric criteria programmed in the NCU 
(ref. 6) and signal validity based on the quality and strength of the received signal 
programmed in the radio receiver. The data recorded on the airplane indicated that 
the Atlantic City (ACY) DME was used throughout most of the test run. However, the 
DME signal was intermittently invalid and caused the Cedar Lake (VCN) DME to be auto- 
matically tuned during the last 49 sec of the functional check. 

Lateral convergence .- The lateral convergence characteristics of the position 
estimate for an offset of approximately 1000 ft for the various position estimate 
algorithm configurations flight tested (i.e., magnitude of Kg, position difference 
solution, and filter time constant) are summarized in table III. The length of time 
of each test and the percent of time that both ILS localizer and DME inputs were used 
to compute the position estimate are also shown in table III. A single offset was 
induced on five runs and two offsets were induced on an additional two runs (denoted 
a and b ) . 

The lateral convergence characteristics obtained from the flight tests were 
similar to those obtained for the same algorithm configuration and filter time con- 
stants used during the simulation tests. The values of the offset obtained at the 
end of each run were generally greater during the flight test than during simulation. 
This was partly because the elapsed time of each flight test was shorter than the 
169 sec used during the simulation test runs. (A tail wind component was encountered 
during the flight tests which resulted in higher ground speeds and subsequently 
shorter times to fly to the runway.) In addition, during the DME signal dropouts, 
the position estimate was computed in the ILX mode. 

Figure 1 1 shows plots of lateral displacement as a function of time for four 
flight test runs. These plots are similar to those obtained during the simulation 
tests. Test run 1 was conducted with the original dual -component position estimate 
algorithm, a 50-sec filter time constant, and Kg set equal to 1. This configura- 
tion resulted in the slowest convergence of the lateral displacement. Test run 2 was 
conducted with the single-component position estimate, a 50-sec filter time constant, 
and Kg set equal to 1 . The higher rate of convergence (compared to test run 1 ) 
shewn on this plot was due to the use of the single-component position estimate 



TABLE III.- LATERAL CONVERGENCE CHARACTERISTICS OF ILS-DME POSITION ESTIMATE 


[TSRV flight test] 


Test 

run 

Velocity 

estimate 

constant 

K 3 

Position 

difference 

solution 

Filter 

time 

constant, 

sec 

Initial 

recovery 

rate, 

fps 

63-percent 

recovery 

time, 

sec 

Offset at 
end of run,* 
ft 

Elapsed 
time of 
test, 
sec 

Time in 
ILD update 
mode, 
percent 

i 

1 

Dual 

component 

50 

18.8 

34.8 

342** 

44* * 

67 

2 

1 

Single 

component 

50 

47.5 

20.5 

-50 

1 22 

97 

3 

1 

Single 

component 

30 

58.8 

14.3 

10 

109 

59 

4a 

1 

Single 

component 

30 

54.4 

14.9 

48 

70 

53 

4b 

1 

Single 

component 

30 

55.6 

14.3 

116** 

33** 

82 

5a 

0 

Single 

component 

30 

61.3 

1 2. 2 

60* * 

46** 

67 

5b 

0 

Single 

component 

30 

56.3 

14.4 

1 05** 

48** 

31 

6 

0 

Single 

component 

30 

58.5 

14.3 

10 

93 

73 

7 

0 

Single 

component 

30 

54.1 

1 5.4 

10 

98 

75 


•Negative offset indicates an overshoot. 

* ‘Insufficient time remaining in test run for lateral convergence to be completed. 


algorithm. However, the lateral error was reduced only to approximately 20 ft and 
then increased and remained at 40 ft until the completion of the test. Test runs 3 
and 6 were conducted with the single-component position estimate and a 30-sec filter 
time constant. For test run 3, Kj was set equal to 1 and for test run 6 it was set 
equal to 0. The lateral convergence was faster on these runs than on run 2 due to 
the use of the quicker time filter constant. The lateral error was nulled smoothly 
(to 10 ft at the end of the test) on run 3. The lateral error on run 6 was also 
smooth but converged to a 35-ft standoff error. This error was maintained for about 
30 sec before converging to 10 ft. This standoff error may be attributed to loss of 
the DME signal, which resulted in the use of the ILX mode. 

The pilot reported no adverse effects to the approach guidance from modifica- 
tions to the position estimate algorithm. The dynamics of the computer -drawn runway 
displayed on the electronic attitude director indicator were smooth with both 30- and 
50-sec filter time constants. 
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Longitudinal convergence.- The longitudinal convergence rates were almost iden- 
tical to those obtained on the lateral axis when the DME signals were valid. When 
the DME signals were invalid, the position estimate was computed with ILS localizer 
signals and the longitudinal recovery rate dropped to zero. 

Figure 12 shows plots of the longitudinal displacement as a function of time for 
two of the flight test runs. Both runs were configured using the single -component 
position estimate with a 30-sec filter time constant. The plateaus where the longi- 
tudinal displacement did not change indicate the times during which the DME signals 
were invalid. In run 3, was set equal to 1 and in run 5a it was set equal 

to 0. It was anticipated that some overshoot would occur in run 3 (Kg = 1 ) because 
of the overshoot obtained in the simulator (table I). However, none of the flight 
test runs resulted in the longitudinal error converging to zero, which indicated that 
a bias error could be present in the DME system. 

In an effort to determine the combined accuracy of the airborne and ground DME 
systems, a comparison was made between the DME value measured on the airplane and a 
calculated distance determined from the location of the airplane (defined by the pre- 
cision radar tracking data) and the known location of the DME ground facility. Slant 
range error due to the difference in altitude between the airplane and the DME sta- 
tion was taken into account. This analysis shewed that the indicated DME values 
measured on the airplane were greater than the theoretical values obtained from both 
DME ground stations used during the test runs. The mean error for all test runs was 
739 ft and the standard deviation was 247 ft. 

The DME ground stations used during these tests were standard radio navigation 
facilities normally used in the national airspace structure. Combined ground and 
airborne accuracy is specified as 0.5 n.mi. or 3 percent of the indicated distance, 
whichever is greater (ref. 7). 

Although the DME system errors resulting during these tests are well within the 
specified accuracy tolerances for en route DME systems, the use of precision DME 
facilities may be required to compute guidance for precision approaches or for super- 
imposition on real-world images (i.e., head-up or head-down displays with electronic 
real-world imagery). The NASA test pilot reported that the longitudinal error was 
only observable on the computer -drawn runway guidance when close to the runway 
threshold (within approximately 1 mile). 


CONCLUDING REMARKS 

The results from simulation tests showed that improvements to the convergence 
characteristics and accuracy of the TSRV area navigation system position estimate 
could be obtained through modifications to the position difference computations, the 
position estimate filter time constant, and the velocity estimate equations. Flight 
test results generally agreed with the simulator results except for the accuracy 
improvements obtained through modifications to the velocity estimate. Accuracy 
improvements with this modification are not evident in the flight data supported by 
radar tracking information. Lateral convergence (influenced by ILS localizer radio 
signals) in nulling an artifically induced position estimate error was substantially 
faster with the single-component position difference solution than with the dual- 
component solution. Convergence in the longitudinal axis (influenced by DME radio 
signals) was approximately the same for both solutions. Use of the 30-sec time con- 
stant filter gains resulted in approximately 66 percent faster convergence than was 
obtained with the 50-sec gains. During simulation tests an overshoot error in the 



position estimate was eliminated by setting an accumulation term in the velocity 
estimate equations equal to zero. 

The guidance drawn on the electronic flight displays was smooth for all position 
estimate algorithm configurations and filter time constants tested. It was noted 
that some longitudinal error could be detected between the computer-drawn runway and 
the superimposed real-world runway displayed on the electronic attitude director 
indicator when the airplane was close to the approach threshold of the runway. Ibis 
longitudinal error was caused by a bias error in the DME indication used in the posi- 
tion difference equations. Although the bias error was within system tolerances for 
en route navigation DME, the use of precision DME may be required for approach to 
landing. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 12, 1 984 



APPENDIX A 


POSITION ESTIMATE FILTER 

The first step in the position estimate process is to develop a velocity esti- 
mate update from the north and east components of the position difference 
(AP n t /AP E t ) obtained from the vector between the last previously estimated position 
and the position defined by the radio navigation information. Thus * 
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where 

K 2 = 900 At/T 2 hr” 1 
K 3 = 1.0 
At = 1 /20 sec 

T = Time constant for position estimate filter 

A velocity estimate is then made by summing the velocity update with ground 
speed obtained from the inertial navigation system: 
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A position update in terms of latitude and longitude is obtained from the system 
velocity and the position difference components as follows: 
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where K =2 At/T. Hiis latitude and longitude position update is based on an 
oblate spheroid Earth model (ref. 3) using the following radii of curvature in the 
meridional (north-south) and normal (east-west) directions: 


( 2 VI 1 

r 1( . = h . + r_ (1 - 2F + 3F sxn (j) . , ~ 

M , t a,t E\ Y e,t-1 ) 6076.1 
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L E \ 
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where 


F = 0.003367 


r E = 20 925 873.79 ft 


The position estimate is found by summing the previous position estimate with the 
position update terms: 
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APPENDIX B 


ILS LOCALIZER-ONLY POSITION DIFFERENCE COMPONENT 


General Solution 


A position difference component may be computed from ILS localizer angular- 
deviation information. Figure B1 shows the geometry of the ILS localizer antenna, 
the runway with an extended centerline, the last estimated position of the airplane 
( <t> e , X ) t _i , and its estimated position on a measured localizer radial. An orthogonal 
coordinate system with its origin at the localizer antenna (<}> , X ) is oriented toward 
true north. ° ° 


A vector DP was found by subtracting the vector \ (defined from the 

t e,t-1 

localizer antenna to the location of the last position estimate) from the vector 

Z (defined from the localizer antenna along the measured ILS localizer radial and 

3T # t — y 

having the same length as vector Z . „ ) . A vector component DP ^ perpendicular 

e,t-1 p,t 

— y 

to the runway centerline was computed from the position difference vector DP^. Then 


DP was broken into north and east components for use in the position estimate 
P/t 

computations . 


Calculation of DP fc 

To determine DP^_ it wa^s necessary to calculate the north and east components 
of position estimate vector Z from the latitudes and longitudes of the last 

6f L“ I 

position estimate and the ILS localizer antenna location: 
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^e,t-1 = (Z e, E, t-1 )i + (Z e,N,t-1 )j 


n.mi . 


* 


A vector z of the estimated position of the airplane along a measured 
localizer radial &as determined in the following manner. Since no radio-based range 
information could be obtained along a radial from the ILS localizer antenna, the 
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magnitude of Z was chosen to be equal to the magnitude of Z . . Ihe length 

r,t e,t-1 

■> -> 

of Z and Z _ was 
r,t e,t-1 


J r,t Z e,t-1 j( Z e,N,t-l) + ( Z e,E,t-l) J 


/2 


n.mi . 


The north and east components of Z were found from the angle d> - ri. 

r, t r t 

->■ 

between the vector Z . and the north axis: 

r, t 

Z „ . = Z cos(d> - tv ) n.mi. 

r,N,t r,t r t 

Z = Z sin(di - tv ) n.mi. 

r,E,t r,t r t 


Thus 


\,t * + <Z r,N,t )j n -" i ' 

Hie position difference vector DP t was found as follows: 

— *• -> -> 

dp. = z - z 

t IT r t t** I 

Calculation of Position Difference AP ,AP 

N,t E, t 

The magnitude of the component DP of DP perpendicular to the runway 

Pit t 

centerline was obtained from the dot product of DP^_ and a unit vector u perpen- 
dicular to the runway centerline. Ihe unit vector u was defined as 


u = -cos ( 4> r )i + sin((J; r )j 


The magnitude of DP was thus 

P*t 


DP = D? • u n.mi. 
p,t t 


= -(Z - Z ) cos ^ + (Z - Z ) sin (|> n.mi. 

r,E f t e,E,t-1 r r,N,t e,N,t-1 r 
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The north and east components of DP were 

P»t 


AP = DP sin <i n.mi. 
N,t p,t r 


AP = -DP cos cb n.mi. 

E,t p,t r 

These position difference components were used directly in the navigation position 
estimate algorithms. • 
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Figure 1 Computer-drawn runway on electronic primary flight display 
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Figure 2.- Geometry pertaining to ILD single-component position estimate: ILS 

localizer antenna, DME, estimated and measured positions of airplane, and runway 
with extended centerline. 




Figure 3.- Components of vector Z triangle used in 

ILD single-component position estimate. 



Figure 4.- Lateral and radial limits of ILS localizer geometry 
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Figure 7.- Plan view of reference coordinate system and initial 

simulation test setup. 











Figure 11.- Lateral component of ILS-DME position estimate. TSRV flight test 
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